A first principles based prediction of electronic and nonlinear optical properties towards cyclopenta thiophene chromophores with benzothiophene acceptor moieties

In the current work, organic cyclopenta-thiophene (CPT) based derivatives (FICR and FICD1–FICD5) were designed by the modulation of end-capped acceptor group of the reference molecule i.e., FICR, to explore their nonlinear optical (NLO) response. The effect of terminal acceptor and donor groups in the tailored compounds was explored by using DFT based quantum calculations. The UV–Vis analysis, frontier molecular orbitals (FMOs), transition density matrix (TDM), natural bond orbitals (NBOs), density of states (DOS), nonlinear optical (NLO) analyses were performed at M06/6-311G(d,p) functional. The LUMO–HOMO band gaps of FICD1–FICD5 were found to be smaller (1.75–1.92 eV) comparative to FICR (1.98 eV). Moreover, the global reactivity parameters (GRPs) were correlated with the results of other analyses. FICD2 and FICD5 with lowest band gap 1.73 and 1.75 eV showed less hardness (0.86 and 0.87 eV, respectively), high softness (0.58 and 0.57 eV−1), and larger absorption spectrum (815 and 813 nm) in gaseous phase and (889 and 880 nm) in solvent phase among all entitled compounds. All the designed chromophores (FICD1–FICD5) demonstrated a significant NLO response as compared to FICR. Particularly, FICD2 and FICD5 exhibited the highest average linear polarizability () [2.86 × 10−22 and 2.88 × 10−22 esu], first hyperpolarizability (βtot) (8.43 × 10−27 and 8.35 × 10−27 esu) and second hyperpolarizability (γtot) (13.20 × 10−32 and 13.0 × 10−32 esu) values as compared to the other derivatives. In nutshell, structural modeling of CPT based chromophores with extended acceptors, can be significantly utilized to achieve potential NLO materials.

emerged as promising candidates for NLO applications, owing to their refractive index dispersion, high structural versatility, swift response rates and optical spectra in comparison to inorganic compounds 7 .In addition, they display great stability in the visible spectrum, natural softness, exceptional optical transparency, flexibility of their chemical structure due to diverse functional groups, wide variety of production processes, electro-optic switching, the generation of terahertz (THz) waves and frequency doubling 8 .Donor-π-acceptor (D-π-A) architecture plays a crucial role in achieving remarkable NLO responses 9 .There is a spatial separation of electrons in D-π-A systems, wherein, rich (D) and electron-deficient (A) moieties are interlinked by a π-conjugated bridge.The π-bridge between the donor and acceptor moieties provides a pathway for delocalized π-electrons 10 .Extending the conjugation length increases the probability of electronic transitions.This architecture exhibits intramolecular charge transfer (ICT) upon photoexcitation, where an electron is transferred from the donor to the acceptor through the π-bridge.This ICT process leads to significant changes in the molecular dipole moment, resulting in large second and third order NLO effects 11 .An effective push-pull mechanism has been established by earlier researchers using of D-π-A compounds tailored with different acceptors 12,13 .These push-pull arrangements affect charge distribution, broaden the wavelength range, intensify asymmetric electron distribution and narrow the energy band gap (E LUMO -E HOMO ), hence increasing the NLO response 14,15 .This architecture is also utilized in the design of organic molecules, polymers, and materials as used in photonics and optoelectronics 16 .Thiophene-based compounds are considered as promising candidates, particularly when engineered with a push-pull architecture 17 .This structural design involves electron-donating and electron-withdrawing groups strategically positioned along the molecular backbone, creating a strong ICT effect.Such a configuration facilitates efficient delocalization of electron density, leading to NLO responses [18][19][20][21] .In current era, non-fullerene (NF) based compounds received a significant attention in the area of photonics due to (i) isotropic electron transference (ii) effective electron delocalization between D and A species, (iii) completely conjugated structure, and (iv) electronegativity 22 .Additionally, NFs materials, characterized by planar structures and tunable energy gaps, demonstrate significantly enhanced stability than other organic systems 23,24 .By utilizing the stronger D and A components in NFs, can enhance the NLO characteristics by developing an appropriate π-conjugated system.
Le and coauthor have designed and synthesized a fused tris-thienothiophene (3TT) building units having molecular packing ability and strong electron-donating effect to achieve efficient NLO properties.The chromophores having three thieno [3,2-b]thiophene (TT) units combined with two cyclopentadienyl (CP) rings, resulting structures exhibit greater planarity, along with an extended electron delocalization length compared to a single thiophene ring.Moreover, increasing the length of conjugation can enhance the NLO characteristics, and a more rigid backbone can prevent bond rotation, thereby improving molecular packing 25 .The structure-property relationship exploited that the choice of ideal length of π-conjugation enhanced the NLO response 26 .Therefore, an elongated π-bridge improves the conjugation, which eventually raises the nonlinearity of the system.Moreover, its light harvesting ability results in red shifted absorption bands 27 .The incorporation of electronegative substituents such as cyano and nitro groups at the acceptor components of NF structures may also significantly influence their electronic characteristics 28 .Therefore, keeping in mind the prior discussion, in this study an organic fullerene free cyclopentathiophene (CPT) based compound is selected as parent compound (FOIC) 29 comprising A-π-A framework.In the parent compound (FOIC), the bulky 1-hexyl-4-methyl groups attached to the π-spacer have been replaced with methyl (-CH 3 ) groups to decrease computational expenses and eliminate steric hindrance arising from the large alkyl chains (Fig. 1).Hence, we designed five compounds (FICD1-FICD5) by altering its architecture from A-π-A to D-π-A and determined their NLO response.The aniline moiety has been selected as a donor owing to its simple, planar, and rigid structure and good electron-donating capabilitywhich enhances the optical nonlinearity of the system 30 .Literature data exploited that this is the first comprehensive DFT work to elaborate the NLO and electronic properties of FOIC based molecules.The main objective of this research

Computational procedure
Current quantum chemical calculations were executed in chloroform media using Gaussian 09 program 31 .For DFT based investigations; molecular structure input and visualizing output files of the entitled systems were generated by using the GaussView 6.0 32 and Chemcraft software 33 .The M06 34 functional along with 6-311 G(d,p) 35 basis set to accomplished current chemical calculations.The M06 functional, coupled with the 6-311G(d,p) basis set, has demonstrated excellent performance and shown a strong correlation with experimental data in numerous studies 36,37 .Therefore, we also used above mentioned functional and basis set in our study.At first, all the geometries were optimized in order to get true minima structures.The NBO 5.0 program 38 was utilized for NBOs analysis i.e., calculating the energy of stabilization, find out intramolecular charge transfer (ICT) among various components (D, π, and A) of reference compound (FICR) and the investigated chromophores (FICD1-FICD5).The TD-DFT 39 is one of the most widely used method for the calculation of excitation energies (UV-Vis) and FMOs analyses because of its reasonable agreement with experimental data.For this reasons, TD-DFT was used for (UV-Vis) and FMOs analyses in this study.The FMOs investigation was performed to calculate the energy gap between HOMOs and LUMOs 40 .The DOS analysis was also carried out to determine the charge density of the studied compounds by using PyMOlyze 2.0 program 41 .All the graphs were developed via using Origin 8.0 software 42 .The other software's such as Multiwfn 43 , GaussSum 44 and Avogadro 45 were employed for the interpretation of results from output files.

Results and discussion
The present research is focused on designing a series of CPT based compounds (FICR and FICD1-FICD5) for NLO materials.For this purpose, already synthesized parent molecule (FOIC) 29 with A-π-A type configuration is taken and by its some structural alteration the reference compound (FICR) is designed as illustrated in Fig. 1.To achieve a simple push-pull, D-π-A architecture, one of the peripheral acceptor of FOIC is exchanged with a donor moiety i.e., 4-methyl-N-phenyl-N-(p-tolyl)aniline, while other terminal acceptor is exchanged as with benzothiophene (BT) based acceptors in order to obtain FICR (Fig. 2).In the designed molecules (FICD1-FICD5), two chloro (-Cl), two nitro (-NO 2 ), two trifluoromethyl (-CF 3 ), two methyl acetate (-COOCH 3 ) and two cyano (-CN) groups are substituted to BT acceptor moiety (Fig. 3).Cartesian coordinates are given in Tables S1-S6.The optimized structures of the entitled compounds are presented in Fig. 4. To study the impact of structural designing on NLO properties, DFT/TD-DFT computations are performed to calculate the transitions of NBOs, spectral absorption, HOMO/LUMO band gaps, ⟨a⟩, β tot and γ tot .It is expected that this research will be significant in the field of NLO and will inspire experimental researchers to synthesize these chromophores with enhanced NLO responses.www.nature.com/scientificreports/

Natural bond orbitals (NBOs) analysis
The electronic charge transfer and intramolecular charge delocalization from HOMOs to LUMOs in the π-conjugated network can be efficiently explained by NBOs investigation 46,47 .It demonstrate efficiency in interpreting orbital interactions within bonds, molecular stability, hyper-conjugative interactions, and hybridization of the studied compounds.Second-order perturbation theory is used to calculate the stabilization energy 48 of the molecules by employing Eq. 1.
Here, i and j indicate donor and acceptor groups, while E (2) is energy of stabilization.Moreover, q i , ε j , ε i and F i,j are donor orbital occupancy, diagonal and off-diagonal NBO Fock matrix elements between the NBOs.Table 1 displays some significant selected transitions while the detailed values of the investigated compounds are illustrated in Tables S7-S12.

Frontier molecular orbitals (FMOs) analysis
An FMO investigation was conducted to compare the electronic properties of the reference molecule (FICR) with its designed derivatives (FICD1-FICD5) along with their GRPs and DOS analyses.It is an outstanding method to determine the electronic transitions, chemical stability and reactivity of compounds 49 .Moreover, the energies of the highest occupied molecular orbitals (HOMOs) and the lowest unoccupied molecular orbitals (LUMOs) were used to determine the GRPs of the studied compounds.The energy gaps of the designed compounds were computed in chloroform utilizing the M06/6-311G(d,p) functional.(Table 2).The other orbitals and their respective energy gaps are illustrated in Tables S13-S18.
Table 2 showed that E HOMOs in FICR and FICD1-FICD5 were found as −5.29, −5.26, −5.33, −5.29, −5.27 and −5.29 eV, respectively.Similarly, calculated values for E LUMOs are as −3.31, −3.36, −3.60,−3.44,−3.36 and −3.55 eV, respectively.In FICD5, HOMO is observed to be similar to FICR but LUMO is at lower energy (−3.55 eV), while with FICD2, both FICR and FICD5 showed HOMO at little bit higher energy (−5.29 eV) as compared to FICD2 (−5.33 eV).The E gap of all the designed compounds (FICD1-FICD5) is found lower than that FICR indicating easily ICT in derivatives.The least energy gap is observed in case of FICD2 (1.73 eV) and FICD5 (1.75 eV) might be due to the presence of strong electron withdrawing moieties (-NO 2 and -CN groups, respectively) on the acceptors as compared to the other designed chromophores.In FICD1, the E gap was observed to be smaller than that of FICR but larger than FICD5 (1.92 > 1.75 eV).The reduction in FICD1 energy gap value as compared to FICR might be due to the replacement of terminal acceptor with benzothiophene acceptor having two -Cl groups.The increase in conjugation as well as inductive effect (-I) and resonating effect (-R) due to chloro group in FICD1 might be the reason of reduction of energy gap in FICD1.The further decreased in E gap in FICD5 might be due to the replacement of -Cl groups with strong electron with drawing -CN groups.In FICD3, the E gap (1.85 eV) was found smaller than FICD4 (1.91 eV) when the -COOCH 3 group was replaced with trifluoro (-CF 3 ) group at the acceptor unit resulting in lesser electron-withdrawing capability.Table 2 showed that, all the designed derivatives have smaller band gap as compared to the reference.The overall descending order of band gap in FICR and FICD1-FICD5 was found as : www.nature.com/scientificreports/For the charge transfer, E gap is a crucial factor i.e., narrower the band gap, larger the charge transfer rate.The FICD2 and FICD5 exhibited the lowest energy gap due to the enhanced conjugation and electron withdrawing nature at acceptor caused by two -NO 2 and two -CN strong electrons withdrawing groups resulting in the most efficient push-pull mechanism out of all the entitled compounds.The data of orbital energies can be used to depict the charge transfer analysis with the help of FMOs pictographs as shown in Fig. 5 and Figure S1. Figure 5 showed that HOMOs were primarily located over donor and π-spacer.In contrast, LUMOs exhibited  www.nature.com/scientificreports/and -CN groups exert strong mesomeric effects, influencing the electronic properties and reactivity of molecules containing these acceptor groups.Thus, a notable transfer of electronic charge is observed from donor to acceptor with the help of π-spacer moiety in the investigated molecules as illustrated in Fig. 5.For HOMO the charge is major located over the donor and π-spacer then in LUMO this charge is significantly moves towards the acceptors.The presence of π-spacer enhances ICT from electron-rich to electron-deficient molecules, demonstrating the effectiveness of these molecular systems as NLO materials.

UV-Visible analysis
In order to estimate the optoelectronic properties of the entitled chromophores (FICR and FICD1-FICD5), their UV-Visible analysis is carried out in chloroform and gaseous phase by employing TD-DFT computations at the above mentioned functional 50 .This analysis provides understanding of charge transference, contributing configurations and types of electronic transitions in the investigated compounds 51,52 .Moreover, the effect of the addition of a donor and modulation of terminal acceptor group on the optical properties of compounds is also studied.By executing TD-DFT computations, six lowest singlet-singlet allowed transitions have been calculated for FICR and FICD1-FICD5 and the results are tabulated in Tables S19-S24.The highest transition wavelengths (λ max ) along with their major contributing percentages are shown in Table 3 and their UV-Vis spectra (gaseous phase) is depicted in Fig. 6.
Table 3 clearly depicts that all the chromophores showed highest absorption in the visible region of electromagnetic spectrum.In gaseous phase, the chromophore FICR exhibited the λ max with major H → L (97%) contribution at the range of 701 nm and oscillator strength (f os ) of 1.77.All the designed chromophores FICD1-FICD5 showed greater λ max values at the range of 734-815 nm than FICR.The compound FICD2 and FICD5 showed the largest absorption wavelength as 815 and 813 nm and the least transition energy value of 1.52 and 1.53 eV with a smallest f os of 0.93 and 0.97, respectively and HOMO-LUMO contribution of H → L (97%).Both AICD2 and FICD5 were observed as the most-red shifted derivatives which might be due to the presence of -NO 2 and -CN groups on the acceptor moiety that reduced the energy gap between molecular orbitals.The λ max (nm) of all the designed molecules is in the following ascending order: FICR > FICD4 > FICD1 > FICD3 > FICD5 > F ICD2.The NLO efficiency of these materials was also calculated from excitation energy.The chromophores with lowest excitation energy exhibit bathochromic shift resulting in efficient NLO properties.The transition energy values (eV) gives the reverse order of λ max : FICR < FICD4 < FICD1 < FICD3 < FICD5 < FICD2.The UV-Visible graph in gaseous phase is displayed in Fig. 6.Table 3. Excitation energies (E), oscillator strength (f os ), wavelength ( ) and contributions of various molecular orbitals of FICR and FICD1-FICD5 in the gaseous phase.MO = molecular orbital, H = HOMO, L = LUMO, (nm) = wavelength.www.nature.com/scientificreports/

Compounds
In solvent phase (chloroform), the chromophore FICR exhibited the λ max with major H → L (94%) contribution at the range of 771 nm and oscillator strength of 2.20.All the designed chromophores FICD1-FICD5 showed greater λ max values at the range of 806-880 nm than FICR (Table 4).Like gas phase, in chloroform media both FICD2 and FICD5 chromophores showed the largest absorption wavelength as 889.29 and 879.95 nm with the least transition energy value of 1.39 and 1.41 eV, respectively and HOMO-LUMO contribution of H → L (94%).Both these compounds have strong electron withdrawing -NO 2 and -CN groups on the acceptor moiety, hence, create a strong push-pull architecture that relatively lower the excitation energy and improved the absorption spectra towards bathochromic shift.The UV-Visible graph in solvent phase is displayed in Fig. 7. Literature study revealed that the polar medium aids in stabilizing the π-π* state alongside the n-π* state by employing a suitable electronic level 53 .Typically, the energy of compound interactions in polar media is governed by non-covalent interactions (NCIs) and polarity influences.This factor delineates dipolar interactions and hydrogen bonding, which play a significant role in stabilizing the molecules' first singlet electronic level.As solvent polarity increases, the molecule experiences a red shift 54 .The excited state is considered to be more polar than the ground state, resulting in greater stabilization of the excited state compared to the ground state in chloroform.Therefore, the current studied compounds showed wider absorption spectra in chloroform than that of gaseous phase which significantly improve the NLO properties of these compounds.

Global reactivity parameters (GRPs)
To explore the reactivity and stability of the designed compounds, the DFT approach is employed to estimate energy values, allowing for the examination of global reactivity parameters such as global hardness (η), global softness (σ), global electrophilicity index (ω), electronegativity (X), ionization potential (IP), electron affinity (EA), and chemical potential (μ) 55 .The GRPs for current study chromophores were calculated by utilizing Koopmans' theorem 56 Equations: S1-S8 given in Table S29.Either donating or accepting capability of a molecule within a system is directly linked to its ionization potential (IP) and electron affinity (EA), respectively.Electronegativity (X) governs a compound's capacity to attract the electron cloud towards its own nucleus 57 .Furthermore, the energy gap is closely linked to the chemical potential, hardness, and stability of a compound, whereas it has an inverse relation with reactivity 58 .
Among all the designed compounds, FICD2 and FICD5 are regarded as the softer, less stable, and most reactive ones owing to their higher global softness value of 0.58 and 0.57 eV, respectively.The global softness of the designed compounds follows a descending order: FICD2(0.58 eV) > FICD5(0.57eV) > FICD3(0.54 eV) > FICD 1 = FICD4(0.52eV) > FICR (0.50 eV) (Table 5).Nevertheless, the trends in global hardness values were noted to

Density of states (DOS)
The DOS analysis was accomplished to validate the findings obtained by FMO analysis for the studied compounds (FICR and FICD1-FICD5) 59 are presented in Fig. 8. On graphs, the HOMOs also called the valence band are represented at right side, whereas the LUMOs known as conduction band are shown at left side.The x-axis depicts the energy gap, which represents the separation of HOMO and LUMO 60 .To explain the findings regarding the density of states (DOS), the compounds are categorized into three segments: donor (donating group), π-spacer (linker), and acceptor (end capped acceptor unit) as demonstrated by blue, green, and red lines, correspondingly.
In the reference compound, donor shows electronic charge distribution of 24.60% toward HOMO and 0.50% to LUMO.Whereas π-spacer shows that electronic charge distribution provides 69.60% to HOMO and 57.50% to LUMO.Additionally, acceptor display electronic charge distribution pattern as 5.70% toward HOMO and 42.00% toward LUMO.On the other hand, in the designed compounds (FICD1-FICD5), donor display the electronic charge distribution pattern as 21.90, 24.

Transition density matrix and exciton binding energy
The TDM investigations plays an essential part in determining the amount of electronic-charge transfer within reference and designed compounds 61 .The TDM is useful for assessing charge density excitation, electron-hole pair localization and delocalization, and excited state interactions between electron-accepting and electrondonating entities 60 .Additionally, it is also helpful in understanding transition behavior from the ground (S0) to excited state (S1) 62 .The TDM of the studied compounds was calculated using the M06/6-311G(d, p) method.The influence of the hydrogen atom is ignored in this analysis owing to its small influence on electronic transitions.The TDM heat map of all the designed compounds demonstrate the nature of transition is illustrated in Fig. 9.In this visual representation, the lower horizontal axis and the left-side vertical axis depict the atom count, while the right-side vertical axis represents the transition coefficient and electron density resulting from interactions with light.Furthermore, the green dots on the dark blue background provide a pictorial illustration of charge transfer.To understand the distribution of charge density, we divided each designed compound into three components: π-spacer (linker), donor and acceptor.As indicated in the pictographs, the highest electronic cloud was apparent on the π -linker in FICR and FICD1 to FICD5.Although a little charge can be observed on the donor region, which successfully facilitate the shifting of electron density towards the acceptor despite trapping by incorporating electron coherence efficiently transferred from the donor to the π-linker.
Moreover, the binding energy (E b ) serves as a crucial parameter aiding researchers in evaluating the optoelectronic characteristics of material 63 .It facilitates researchers in computing the coulombic interaction forces between the hole (h) and electron (e).At excited state, exciton dissociation has an inverse relationship with the E b , however there is a direct correlation between the E b and the coulombic interaction between the h and e. Reduced E b results in higher charge mobilities 64 .The E b of designed chromophores FICR and FICD1-FICD5 was determined using Eq. 2 65 .
(2) www.nature.com/scientificreports/In Eq. ( 4), E H−L signifies energy difference between HOMO/LUMO.E opt represents that the minimum energy required for the initial excitation from S 0 to S 1 , resulting in generation of electron and hole 66,67 .Table 6 displays DFT computed findings for E b of the above-mentioned chromophores.
The E b value revealed that, FICD2 and FICD5 demonstrated a greater extent of charge separation when compared to other derivatives.Table 6 presents the calculated values of the HOMO-LUMO gap (E H -E L ), the first single excitation energies (E opt ), and the excitation binding energies (E b ), combined with its high current charge density, makes FICD2 and FICD5 a highly promising candidate for achieving exceptional efficiency in NLO compounds.The observed correlation between binding energy and TDM outcomes demonstrates a remarkable degree of charge transfer within molecules.The following compounds' binding energy values are listed in decreasing order: FICD4 > FICD1 > FICR > FICD3 > FICD2 > FICD5.www.nature.com/scientificreports/

Nonlinear optical (NLO) analysis
Linear polarizability (α) defines the potential of electric field to distort the molecular electron distribution 68 .The hyperpolarizability (β tot , γ tot , etc.) is a key parameter that characterizes the nonlinear optical behavior of molecules and atoms, encompassing a wide range of complex phenomena within the field of nonlinear optics.The first hyperpolarizability (β tot ) is ascertained through the finite field method, where the proximity of an applied electric field influences the system's molecular energy.The relationship between molecular structure and NLO effects is defined by the charge transfer capability and hyperpolarizability of a substance.When a material displays a substantial NLO response, it demonstrates an elevated dipole moment, significant linear polarization, and substantial hyperpolarizability values.The dipole moment (µ tot ) 69 of all the entitled compounds was determined by using Eq. 2.
Similarly, the linear polarizability ⟨a⟩ 70 was computed with the help of Eq. 3.
The magnitude of first hyperpolarizability β tot 71 was measured with the help of Eq. 4.
where The second hyperpolarizability γ tot 72 was calculated by using Eq. 5.
The calculated results for the electric dipole moment (µ tot ), linear polarizability (< α >), total first order hyperpolarizability (β tot ), and total second-order hyperpolarizability (γ tot ), along with their respective tensors, have been presented in Tables S25-S28 for the FICR and FICD1-FICD5 systems.A summary of the key discoveries can be found in Table 7.
The dipole moment arises because of variations in electronegativity.A larger disparity in electronegativity corresponds to an increased dipole moment (µ tot ).In the context of the FICD5 molecule, it exhibited the highest total dipole moment (µ tot ) of 19.76 Debye (D) because of the notably high electronegativity cyano groups among all derivatives.The FICR and FICD1-FICD5 have shown the trend for dipole moment values as; FICD1 < F ICD3 < FICD4 < FICR < FICD2 < FICD5.The increased electronegativity of atoms results in the attraction of charge density towards them, thereby leading to increase of linear polarizability.Consequently, the greatest linear polarizability, measured at 2.88 × 10 −22 esu, was investigated in FICD5.This phenomenon can be associated with the incorporation of electronegative functional groups within the acceptor moiety.The observed linear polarizability trend of FICR and FICD1-FICD5 is FICR < FICD1 < FICD3 < FICD4 < FICD2 < FICD5.In the context of the FICD5 model, the dominant tensor component responsible for linear polarizability was oriented along the x-axis, denoted as α xx , with a magnitude of 5.94 × 10 −22 esu.This specific tensor component plays a prominent role in determining the overall linear polarizability value.
The analysis revealed that transitioning from the A-π-A configuration to the D-π-A configuration, coupled with an increased electron-withdrawing effect, induced pronounced charge transfer and a redshift in the spectral absorption, leading to enhanced optical nonlinearity in the derivatives.Notably, the narrowing of the band gap resulting in enhanced NLO property.Both the band gap and NLO property having inverse relation where compounds with narrower band gaps typically demonstrate superior NLO performance.It is noted that, both the compounds FICD2 and FICD5 having the lowest energy band gaps (1.73 and 1.75 eV respectively) and highest 2nd hyperpolarizability values (13.2 × 10 −32 and 13.0 × 10 −32 esu respectively).Due to their substantial polarizability and hyperpolarizability values, these designed compounds hold promise as prospective candidates for NLO materials.Among them, FICD2 and FICD5 emerge as a particularly outstanding constituent for nonlinear optics.

Conclusion
Herein, a series of CPT based derivatives (FICD1-FICD5) was designed from a NF molecule (FICR) by modifying the configuration from A-π-A to D-π-A.It is observed that acceptors pay a notable impact on electronic properties and to enhance the NLO response of the tailored compounds than the reference compound.Our performed analyses demonstrated that all proposed molecules experienced efficient charge transfer from donor to acceptor across a π-bridge, resulting in a lower energy gap (1.75-1.92eV) than FICR (1.98 eV).The TDM and DOS data also supported the charge transfer process between HOMO and LUMO with an effective way.The FICD1-FICD5 showed greater absorption and less transitional energy values as compared to FICR.The broader red-shift spectrum value (λ max = 815 and 813 nm in gaseous phase and 889 and 880 nm in solvent phase) was exhibited by FICD2 and FICD5 as compared to all the other designed compounds.Moreover, NBOs analyses revealed that prolonged conjugation causes stability in the studied molecules.The significant NLO properties in terms of < α > (2.86 × 10 −22 and 2.88 × 10 −22 esu for FICD2 and FICD5, respectively), β tot (8.43 × 10 −27 and 8.35 × 10 −27 esu for FICD2 and FICD5, respectively) and γ tot (13.20 × 10 −32 and 13.00 × 10 −32 esu) were investigated as compared to FICR and other designed compounds.The promising NLO results for FICD2 and FICD5 might be obtained due to the strong electron-withdrawing groups.Effective NLO responses of current study chromophores may serve as a driving force for researchers to synthesis such efficient NLO chromophores for photonic applications.
FICD2 (0.86 eV).The ΔN max values reveal that the FICD2 and FICD5 compounds have the highest value at 5.18 and 5.06 eV, demonstrating a constant charge transfer throughout these molecules.The FICD2 and FICD5 were predicted to be the most favorable compound as compared to the others owing to their maximum softness value of 0.58 and 0.57 eV and lowest hardness value of 0.86 and 0.87 eV, respectively.Compounds (FICD2 and FICD5) with lower softness values exhibit higher polarizability with grater ICT might be consider as efficient materials for NLO applications.